On the Interpretation Of the Local Dark Matter 
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ABSTRACT 

The cause of the extended rotation curves of galaxies is investigated. It is shown that con- 
ventional sources and most exotic sources for the needed gravitational fields are implausible. We 
suggest spatial fluctuations in a scalar field, similar to the inflation field, are responsible for the 
gravitational fields. These fluctuations play the role of 'dark' halos around galaxies. They take 
~ 10^ yrs to develop and could not have been important in the early days of the universe. When 
galaxies are clustered, a A— term appears naturally in this theory. The universe's present energy 
density associated with these scalar field variations isf2~l/2 — 2/3. A possible scenario is 
suggested in which the cosmic scale factor TZ{t) would have experienced a recent acceleration. A 
discussion of further observations and theoretical work needed to resolve some ambiguities in the 
theory is given. 

Subject headings: galactic rotation, dark matter, inflation, dark halo, quintessence 



1. Introduction 

This paper speculates that the local dark mat- 
ter found in halos of galaxies may actually be 
fluctuations in a scalar field. The general distri- 
bution of dark matter has been summarized by 
Trimble (1985); Ostriker (1993). By 'local dark 
matter' (LDM) we refer to the unexplained source 
of the extended rotation curves, Vrot{r) vs. r, of 
the Milky Way galaxy (Fich & Tremaine 1991) 
and other galaxies. This definition means that the 
LDM may not be the sole component of the dark 
matter considered in cosmology. But we deal with 
this restricted class because the Sun is on the 'flat' 
portion of the Milky Way's rotation curve and we 
can calculate the LDM's density in our vicinity; 
from this one can derive severe constraints upon 
the halo's possible composition. 

The rotation curves of dwarf, elliptical, and spi- 
ral galaxies have been reviewed by Casertano & 
van Gorkom (1991), de Zeeuw & Franx (1991), 
and Sofue & Rubin (2001). While a correlation is 
seen between the shapes of the rotation curves and 
the surface brightness and morphology of galax- 
ies, modeling them requires a mixture of baryonic 



matter and LDM for which the effective density 
falls off as r~^ at large r. 

So, rotation curves are generally observed to be 
composite with the LDM contribution dominating 
the rotation curves of normal spiral galaxies out- 
side their central core regions; see Rubin (1980, 
1982, 1985). They may be schematized by a lin- 
ear rise of the rotation velocity V to a, value Vh 
at a galactocentric distance r = Rh, followed by 
a constant rotation velocity (± ~ 10 km s~^) for 
Rh < r < Ri; here Ri is set only by observa- 
tional limits of finding radiant baryonic material. 
About 74% of field spirals show these simple rota- 
tion curves (Sofue & Rubin 2001). 

We will ignore those galaxies in which large 
baryonic contributions to the total mass make it 
difficult to extract the LDM's effect on the ro- 
tation curve. These include: galaxies with pe- 
culiar shaped rotation curves in r < Sofue 
et al (1998); compact bright galaxies which may 
have declining velocity curves for r 2 — 3Rh 
Casertano & van Gorkom (1991); and low sur- 
face brightness galaxies for which Vh is small < 
100 km s~^, and there is clear evidence of large 
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amounts of baryonic material at r > 20 kpc, (Pick- 
ering et al 1997; Quillen & Pickering 1997; Mc- 
Gaugh, Rubin, & Blok 2001); we regard these as 
nascent galaxies. 

Observations of the rotation curve of the Milky 
Way galaxy (Clemens 1985; Merrifield 1992) show 
it is composite, with this typical underlying LDM 
distribution. For our galaxy, adopting a solar 
galactocentric distance of Rq ~ 7.8 kpc (Fich 
& Tremaine 1991), one has Rt — 5 kpc and 
Vh ^ 200 - 210 km s"^ The hmiting value Ri for 
which y = 14 is uncertain; it is generally assumed 
to be at least 25 kpc and is more likely to be at 
least 50-60 kpc and may well be >~ 200 kpc (Fich 
& Tremaine 1991; Kulessa & Lynden-Bell 1992). 
In some spirals, one actually observes Ri > 60 kpc. 
If the LDM is the same as the dark halo material 
used in galaxy-galaxy gravitational lens models, 
then Ri > 200 kpc is appropriate (Guzik & Seljak 
2002). We shall adopt Ri > lORh in calculations. 

For other spiral galaxies, generally Rh has 
about the same value. However while Vh often 
is ~ 200 km s""'^, values in the range ^ 100 to 
250 km s"^ do occur. There also is a smaller class 
of galaxies of interest in which the rise is more 
gradual, Rh ~ 10 kpc, with the rise from r = 
to r = Rh better represented by a convex curve, 
see examples in Persic & Salucci (1995); Generally 
these have not been observed much beyond Rh- 

The Milky Way schematized velocity curve, 
assuming spherical symmetry, gives for r < 
Rh, (p) = 6.4 X 10-24 g cm-3 = 0.086 M© pc'^ = 
4.0 hydrogen atoms cm""^. For r > Rh one has 
Pldm = \{p){Rh/''')'^ ■ At the Sun's galactocentric 
distance, i?o, one has po = 0.87 x 10~2'*g cm~^ = 
0.012 Mq pc~^ = 0.54 hydrogen atoms cm~^. 
Masses of interest are: M{Rh) = 4.5 x 1O^°M0; 
M{r > Rh) = {r/Rh)M{Rh); and M(i?o) = 
7.0 X lO^^M© = Mg- For r < Rh the calculated 
mean density may not be too useful for represent- 
ing the LDM because only a contribution for the 
central spheroid bulge has been subtracted and 
an uncertain stellar disk contribution has been ig- 
nored. In this region, however, the density of the 
LDM must fall appreciably below the r"^ law ( 
for otherwise the horizontal portion of the veloc- 
ity curve would continue inward.) The value of 
Rh depends somewhat upon this uncertain bary- 
onic disk contribution to the rotation curve for 
r < Ro (Kuijken & Gilmore 1989a;b;c). However, 



for r > Rq, the calculated masses and densities 
are dominated by the LDM contributions ( e.g. 
Merrifield (1992)) and using these calculated den- 
sities and masses for representing the LDM cannot 
overestimate the true values by more than a factor 
~ 2 if the rotation velocity is to have the observed 
value and lack of r— dependence. Similarly the 
value of Rh may be somewhat underestimated. 

In our vicinity where p oc r~^, the LDM col- 
umn height, is /q°° n dz = 7ri?oWo/2 = 2.0 x 
10^2 hydrogen atoms cm"^ = 145M0 pc"^; up to 
1 kpc the projected surface density is 11.8 Mq pc"^ 
If the LDM were hot gas, the emission measure 
would be J^n^dz = -kRotiq/^ = 1.8 x 10^ cm~^ 
pc. 

In the next section, we show that no population 
of observed astronomical material can play the role 
of the LDM. We also show that proposed exotic 
sources, massive neutrinos, primordial mini-black 
holes, black fluids, etc. are highly unlikely to be 
the LDM. 

In the third section, a scalar field source term 
is examined which produces gravitational acceler- 
ations similar to that observed. The fourth sec- 
tion contains proposals for observational tests and 
a discussion of some uncertainties in the theory 
that still need to be resolved. 

2. Unacceptable LDM Sources 
2.1. No Interstellar Constituants 

The observed LDM mass density for r = Rq is 

very much higher than the observed energy densi- 
ties of magnetic field, cosmic rays, and radiant en- 
ergy; these cannot be primary constituents of the 
LDM. The LDM cannot have a significant compo- 
nent of ordinary interstellar gas, for the column 
densities needed are much higher than those ob- 
served; see pg 525 of Allen's Astrophysical Quan- 
tities (Cox 1999), hereafter referred to a.s AQ. A 
very high temperature medium is excluded by the 
low emission measures observed in the x-ray region 
(Marshall & Clark 1984; McCammon & Sanders 
1990). One cannot make the LDM out of, say 
^5 — 10 Mg of He gas, without providing an ex- 
planation as to where the missing ^ 15 — 30 Mg of 
H, expected by observed abundances, disappeared. 
Nor can the LDM be explained by large dust parti- 
cles. By weight, dust is comprised of 'metals' with 



2 



universal abundance ( by weight) Z < ^ 0.02. For 
every 1 Mq of LDM attributed to dust one must 
account for a missing 50 Mq of H & He. Con- 
sequently all particulate matter up to the size of 
Uranus and Neptune can be excluded, as well as 
very small compact molecular clouds (Alves, Lada, 
& Lada 2001). 

2.2. No Stellar Constituants 

The mass distribution in the solar vicinity is 

given by Table 19.9 of AQ. The only low lumi- 
nosity stars that could contribute significantly to 
the LDM are some low mass main sequence stars 
and white dwarfs/neutron stars. Looking at edge- 
on galaxies, the surface brightness for any class 
of stars with a constant L/M should fall off as 
1/r: this is not observed implying only very low 
luminosity stars could be a main LDM eonstitu- 
ant. Trimble {AQ, p. 530) gives as an lower limit 
for our galaxy, M/L ^ 18 to r = 35 kpc; this cor- 
responds to stars with M~ O.IM© {AQ, p.487). 
But the integrated stellar mass function up to 
M= 0.1 - 0.2 Mq is only 22% - 35% of the to- 
tal main sequence stellar mass function {AQ, pg 
488). Therefore, if the LDM of ~ 4 Afc seen up to 
r = 35 kpc is due to such low mass stars, one must 
account for a missing 12 — 20 Mq of main sequence 
stars associated with them by current ideas of star 
formation. Similarly, white dwarfs and neutron 
stars are excluded by the same type of argument. 
In their formation at least a comparable mass of 
material is expelled into the interstellar medium 
and it is not observed. 

Now, these arguments cannot exclude the pos- 
sibility that a very old population of stars was 
formed (or captured) with an initial mass function 
favoring objects of 1 — 100 MjupUer or that popula- 
tions of small black holes were somehow produced 
ab initio. Normally, such proposals are not de- 
tailed enough to discuss either the lack of expected 
associated matter or the type of physics needed to 
constrain the mass range. It is also possible that 
massive mesons or some other exotic particles are 
present. These possibilities can be restricted by 
the more general physical arguments, given in the 
next subsections. 



2.3. 'Dark Fluid' Hydrodynamics 

Suppose; the; c;nergy-momentum tensor of the 
LDM can be represented by a simple fluid with 
scalar pressure, P = P{p)- First, assume hy- 
drostatic equilibrium. Then, for the outer region 
r > Rh, using p oc r~^, one can solve for an equa- 
tion of state. One finds P = \V^p + Pq where 
Pq oc p^h.l'^'? ^ constant , so that for r > Rh we 
have an isothermal sphere solution. For finite p 
at r = 0, the isothermal Bonnor-Ebert solutions, 
Ebert (1955); Bonnor (1956), discussed by Alves, 
Lada, & Lada (2001), are appropriate . They have 
too small a range in which p oc r~^. Also, for our 
thermal velocity, ~ Vh, they are Jeans unstable 
for i?/, >~ 20 kpc and hence not acceptable. So, 
hydrostatic equilibrium with a simple equation of 
state, is unlikely. But, this argument would not 
necessarily exclude a gas in a state of partial ion- 
ization. 

Non-hydrostatic equilibrium models can be re- 
jected. Postulating centrifugally supported ellip- 
soidal LDM matter runs into another serious prob- 
lem. Suppose V^/r = (1 — ^)V^/r for sufficiently 
large r > Rh; then the effective thermal velocity 
is reduced by a factor However, one would 
then have to explain why the total LDM angular 
momentum increases cx r^. 

Turbulent support by means of a term v • 
Vvr is unlikely. Mass conservation requires 
4:'KRlp{vr) <~ 1Mg/10^" yrs; this gives as an 
upper limit (vr) ~ 1 km s~^. It would be difiicult 
to maintain this low average radial velocity over 
the entire circumference at r = i?o unless (w^) is 
small, say~ (5 km s~^)^. Using (v^)^^^ as rep- 
resentative turbulent velocity, one would require 
a characteristic turbulent length scale <^ 5 pc 
to balance the gravitational acceleration / Rq. 
Again, a mechanism for providing this small scale 
would be required. 

Since p oc is observed, mass conservation 
would limit outward (or inward) streaming of mas- 
sive particles to (vr) ~ 1 km/s~^. This is an un- 
likely restriction on any hypothesized relativistic 
particles. 

2.4. Galaxy-Galaxy Collisions 

Galaxy-galaxy collisions involving halos pose 
difficulties for models of the LDM involving parti- 
cles not interacting with ordinary matter nor ca- 
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pablc of radiating. If the particles undergo col- 
lisions among themselves, then the hales should 
relax adiabatically after a galaxy-galaxy collision, 
not isothcrmally (as is observed). If particle- 
particle collisions are unimportant (as e.g. mini- 
black holes), then the particles' average velocity 
must be of order Vh or less for them to remain 
bound to a galaxy. Then if two galaxies with halos 
collide, one knows that the particles' trajectories 
would be modified by the non-sphcrically symmet- 
ric gravitational potential present during the long 
galaxy collision. After the galaxies separate, there 
would be no restoring force to re-establish spheric- 
ity for the halos, and the time evolution of the 
particles' distribution function would be governed 
by Liouville's equation. 

This is a good argument which is not yet conclu- 
sive because we do not have observational studies 
of halo asymmetries at large r. In principle it is 
possible to study the aftermaths of halo collisions. 
About 5% of the galaxies are in rich clusters of 
galaxies (see Bahcall AQ, p. 613) and in the core 
regions , --^ 1.5/i_i Mpc radius, the average galaxy 
halo with Ri > 50 kpc experiences at least one 
halo-halo collision in 10^" yrs with a short colli- 
sion time < 10^ yrs, the collision frequency being 
oc Rf. Similarly, In galaxy groups and poor clus- 
ters, which contain ~ 50% of the galaxies, the av- 
erage halo with radius Ri > 250 kpc experiences 
collisions in 10^" yrs. In some groupings, of course 
such as the Milky Way-M31 Local System even 
smaller halos should show severe distortion if com- 
posed of particles. Field galaxies should generally 
keep undistorted their primitive halo structure. 

Assiiming halos show approximate spherical 
symmetry out to ^ 50 kpc , we need conclude that 
any exotic LDM candidate must be able to radiate 
efficiently during galaxy-galaxy collisions, either 
at presently undetectable wavelengths {i/ < 10 
Mhz) or in very low energy particle emission ( e.g. 
pairs of neutrinos) to avoid adiabatic relaxation 
in collisions, and that, if considered as a fluid, the 
equation of state cannot be reduced to a simple 
P — p relation if gravitational instability is to be 
avoided. 

3. The Scalar Field 

As an alternate to specifying radiative prop- 
erties and thermodynamics of unknown parti- 



cles, we suggest that a class of fluctuations in 
a scalar field (p form gravitational potential wells 
into which baryonic matter may flow, forming lu- 
minous galaxies in their center regions. These 
potential wells have gravitational mass and are 
the dark halos. There is no reason to assume 
all halos have luminous galaxies associated with 
them, or that all luminous galaxies are associated 
with these dark halo potential wells. 

The halo potential wells will be specified by 
the energy momentum tensor Tab determined by 
= -V{-g){Tab/2)Sg^'' once a La- 
grangian jC{(f>) is chosen. The associated field 
equation is given by 5C/5(j) = 0. We use as 
a Lagrangian C = —jO^L'^+K where = 
(I>,i4',j + mifnjlcjy^il — ^(p^) — A]. Here rrij is an 
assigned timelike vector, m°'ma = rm?, and the 
coupling constant is the dimensionless a^, not the 
usual n = SttGc^^ used for fluids. If one sub- 
stitutes \I' = a0, and A = A = 0, then £ is 
the Ginzberg & Landau (1950) Lagrangian used 
in studying phase transitions and the Higgs field 
and introduced into astronomy for the inflation 
field by Guth (1981). The parameters A & A do 
not affect the behavior of (f>; they are introduced 
to allow us to shape the extent of each individual 
dark matter halo. The energy-momentum tensor 
then is Ty = a^(Ly — ^gijLl) + QijA. The form 
of the wave equation in fiat space is: 

a,V-V2.^-m2</,(l-<A2)=0. (1) 

Since T^t oc (fy^, the requirement that a field the- 
ory has an effective mass density falling off as r^^ 
at large r and finite at the origin, basically forces 
the use of form of the wave equation in fiat space 
to be that shown, restricting permissible forms of 
C. The wave operator gives the behavior for 
the effective density at large distances; the sign of 
the 'mass' term m."'m,a is opposite that used in the 
Klein-Gordon equation because we need solutions 
finite at the origin; and a non-linear term must be 
introduced in the potential term to limit growth. 
We use a scalar field because it is simple and it is 
conceptually economical to see if a descendent of 
the infiation field can play a contemporary role. 

3.1. Approximate Halo Solutions 

We summarize the discussion of the fiat space 
solutions given in the Appendix. We assume that 
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the field always has a background of infinitesimal 
fluctuations, similar to those seen in the CMB, and 
consider only solutions which could have grown 
from these fluctuations 0(r, t = o) ^ and re- 
main finite. The nonlinear cubic term has two 
significant functions: (1) it limits the growth of 
timc-dcpcndcnt solutions; and (2), it may provide 
a strong interaction between two sets of fluctua- 
tions. Aside from these, the term has only a mod- 
est effect on the analytic form of the two classes 
of solutions found. 

The first class of solutions are the large 
wavenumber modes, > rri^. They represent 
traveling waves, t-waves, for short, whose frequen- 
cies are amplitude dependent. [ Wave packets of 
the t-waves are problematical; they would act as 
tachyons since the group velocity du/dk = k/co is 
greater than that of light.] These t-waves are sta- 
ble, non-localizable, and remain infinitesimal for 
f > 0; for this reason they will be mostly ignored 
in this paper. 

The second class of solutions, composed of 

modes with fc^ < m^, do not 'travel' and can grow 
to finite amplitude. With appropriate boundry 
conditions, they exhibit anharmonic periodic mo- 
tion, first growing exponentially fast from a small 
initial value until (p"^ = {(p"^) < 1 ; then de- 
acceleration takes place and the solutions reach 
their maximum amplitudes and then decrease once 
more to their initial amplitudes. The rise time 
from the fluctuation level seen in the cosmic mi- 
crowave background, to the maximum value is 
~ 1 — 5 X 10^ yrs for ~ ^ — 1, and ~ 
3 kpc ^ 10'' yrs (see below). They could not have 
arisen early in the history of the universe. They 
are standing waves (or s-waves.) 

We suggest these represent galactic halos. Be- 
cause their time derivatives make relatively small 
contributions to the energy-momentum tensor, 
(see section A. 3) time variations for the s-waves 
can usually be neglected and the steady state so- 
lution, k^ m is our model for a typical galactic 
halo. 

The approximate steady state solution ^ is, us- 

^We use rrir = mg = = and t for ct in formu- 

lae. We regard as small, and ignore the small {dtcf))^ 
terms and terms of order ~ VJ^/c*. The time depen- 
dent s-wave halo solutions are of the approximate form: 
<j)=[l + /i2sine(/i2,t,r)]5 05 where Q<h?{r) <1 and 



ing rh = m^y{l — a^) and equation A6: 

Sm TTIT 

dig ~ a — ; ~ a for mr < 7r/2; (2) 

mr 

sinm(r — a^7r/4) sin(mr) „ 
cpa =i a K, — ^ for mr > w. 

mr mr 

(3) 

Here, a is an arbitrary amplitude < 1. 

One may use the steady state solution as rep- 
resentative for evaluating the energy momentum 
tensor for s-waves, using the approximations of 
equations (2) & (3). For the radially symmetric 
interior Schwarzschild metric, dr^ = B{r)dt^ — 
A{r)dr'^ — r'^dVi^, standard relations ( see, e. g. 
Weinberg (1972)) give: 

rA~'^ = f-j Mo + A]r^ dr 

\nAB = J [2Kpa]r dr + constant; 
f = -Tl, = -{2A)-'drB. (4) 
where 

Kpo = ^a'[(mV + (9r'A)')-Am'], (5) 

and c^f is the radial gravitation acceleration a 
non-relativistic particle experiences because of this 
s-wave solution. 

For the case A = A = 0, 

r(l - A"^) = ^a^a^r[l - sin^ mr/{mrf] (6) 

/ = -a''a\2r)-^[l-sm2mr/{2mr)] = -{G/^)Mhaiolr'' . 

(7) 

For spherically symetric metrics, there are two dif- 
ferent definitions of the mass if AB ^ constant; 
the first is the volume integral of the energy den- 
sity ( found in the solution for A) and the second 
is the effective mass determining the differential 
acceleration of a test body, / (found in the solu- 
tion for dr-B). We give mass formulae according 
to the second definition. Outside a cluster, where 
AB = 1, the mass defintions are equivalent. 

The parameters of the theory are easily deter- 
mined. The rotation velocity, — {V (r) / c)"^ / r is 
determined by /. Since Rh — 5 kpcs for the 

is an elliptic integral (See the Appendix). 
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Milky Way, one determines =3.2 pes by set- 
ting 2mRh = TT . Similarly, one finds a^^a^ = 
2(V?i/c)^ = 8.8 X 10"'^; from the discussion in Sec- 
tion 3.2, we suggest ~ 5 ^ 3- This choice is 
also suggested by the fact that the maximum ob- 
served rotation velocities of galaxies, correspond- 
ing to a ~ 1 is less than twice that of the Milky 
Way. If the stellar and LDM contributions to a 
galaxy's rotation curve can be separated in the 
region r < Rh, one could separately determine 
& from d Vrot/dr and in cases where 
1, since the acceleration depends upon m 
rather than on m there. 

3.1.1. Isolated Physical Halos 

The infinite halo solution shown above is for 

one halo in an otherwise empty universe. It needs 
modification because exceedingly small values of 
(j) should not significantly contribute to the energy 
momentum tensor. There is no obvious physical 
interpretation for such contributions when they 
produce densities much less than mean baryonic 
densities in the neighboring intergalactic medium, 
or much less than the critical cosmological density 
Pc- In studying halo stuctures, we may supress 
these low background density contributions by in- 
troducing a reference energy level A = Aq > 0, 
or a related limiting value r = Ri, by impos- 
ing other physics. One introduces a 'cut-off' for 
the halo Tat by chosing A to require that at 
r = Ri = (a/m) Aq^^^ one has Kpo = and 
dln(AB) /dr\r=Ri = 0. Then, we may join an ex- 
terior Schwarzschild solution for r > i?; to the 
interior halo solution, in effect regarding Ri as the 
end of the physical halo. Suppose mr ^ 1. The 
halo mass interior to Ri is 

{G/c^)M{Ri) = {l/6)a^a^Ri, (8) 

The gravitational acceleration for the interior re- 
gion r < Ri is then given by 

/ = -a'^a'^{2r)~'^[l-sm2mr/{2mr)]+rXoa^m'^/3. 

(9) 

So that / ^ -{G/(?)M{Ri)/Rf ,as r ^ Ri; for 
r > i?;, one has / ~ —a^a^Ri/Qr^. This rep- 
resents an isolated physical halo appropriate for 
field galaxies. Since the quantity / is an ob- 
servable, one should see that the rotational ve- 
locity does decrease in the outer parts of halos, 



V{r) = Vh[l - f (r/i?i)2]i/2 for r < Ri; again ob- 
servations could determine i?;. [Note that we are 
explicitly ignoring contributions to / from 

any more distant 'point' sources.] 

For this case, when a halo is limited by the den- 
sity of surrounding baryonic matter, a very crude 
estimate for i?,; can be made. Suppose we con- 
sider a halo galaxy surrounded by baryonic dwarf 
galaxies, each of mass 10® — 10^ Mq with a mean 
spacing s = 120 — 250 kpc. [ We used fldwarf = 
^baryon/3 ~ 0.01 Bristow & Phillips (1994), a 
cosmological critical density pc = 2.0h^ x 10~^®, 
and = 0.5 in making these estimates.] For 
the Milky Way-like galaxies the halo density be- 
comes comparable to the dwarf galaxies' mean 
baryonic density when Ri ^ s. Since dwarf galax- 
ies tend to cluster around large galaxies, this seems 
to be a reasonable order-of-magnitude guess for 
Rl when intergalactic baryonic matter is present. 
This gives M{s) = 5—10 Mq, values which are 
1/3 of that given by equation (7) (when A = 0). 
Otherwise, an upper limit Ri ~ 1.6 Mpc for 
Milky Way-like field galaxies results when the halo 
density becomes comparable to pc- The uncer- 
tainty in specifying A (or Ri) does not affect the 
value of / in the inner parts of halos. 

3.1.2. Clustered Physical Halos 

These formulae hold only when other halos are 
not close-by, i.e. Ri <C .Ro where Rq is half the 
mean spacing between halos (galaxies); we refer to 
this situation as case I solutions. For field galax- 
ies Ro ~ 2.2 Mps. In crowded regions such as 
clusters, one has another situation (case II solu- 
tions) where Ri > Rq: the halo boundry must be 
re-examined because the 0— field may make sig- 
nificant contributions to the energy density there. 
We find it necessary to introduce a A 7^ term in 
order to define reference background levels appre- 
ciably greater than zero. 

In the next section we note that a solution for 
(j> can be broken up into many effectively inde- 
pendent parts if their centers are far from one 
another. Assuming this, one can consider that 
(j) can be represented by cluster of similar dark 
halos with mean spacing between halo centers 
r = 2i?o ^ m^^ . Represent each halo solution by 
a cell with a central 'bump' on top of a ' plateau', 
the plateau of one cell joining smoothly onto the 
plateaux of the adjacent cells. In the weak gravi- 
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tational approximation wc can treat each cell sep- 
arately. Consider two adjacent halos with cen- 
tral amplitudes oi, 02- Because we are limited by 
our spherical symmetry assumption, look only at 
the forces and densities along the lines connect- 
ing their centers. Then, at the join point be- 
tween cells, r = ri, the individual halo densities 
Kphaio = («^/2)[ai?']~^ ~ Am^] + A match, provid- 
ing af j-j^^ = a2r^^, where ri + r2 = 2/?o ( and 
mri, mr2 3> 1). Wo can then regard this location 
as defining the edges of the halo bump's interior 
Schwarzschild solution, by forcing dlnAB/dr = 
there with the choice Xrn^ = a^r\'^: an exte- 
rior solution holds for r > vi. The underlying 
background density (our 'plateau') will be repre- 
sented by the; A term; because of the local spher- 
ical symmetry assumption A is constant in a cell. 
If now we choose as the uniform background den- 
sity A = o?\rr? jl we can restore a halo's interior 
density phaio to be that of the 0— field, its physical 
value. One has: 

/ = -a2[a2/2r-AmV3]+Ar/3-*0 (10) 

as r — > ri. For r\ = r2 = Rq one finds each halo 
bump has a mass {G / c^)M{Rq) = [l/^ja^a^RQ, 
superimposed on the mean background density 
level PA = \{phaio) = (3/87r)M(i?o)/i?g. 

In each cell Ar/3 acts like a differential tidal 
force; the A— term causes the entire cluster to ex- 
perience an expansion force. Our choice for A re- 
quires that the total acceleration on a test mass 
vanish at cell interfaces. 

For example, suppose we have 125 galaxies like 
the Milky Way, with {a?)a'^ = 8.8 x 10"'^, in a clus- 
ter with radius Rd = 1.5 Mpc. [ We ignore the 
possibility of dwarf baryonic galaxies in the clus- 
ter limiting the value of Ri] Then Rq = 300/i~^ 
kpc. The typical galaxy halo mass is ^ 13ft-~^ Mq 
and the total mass of the cluster dark matter is 
~ 1.7 X 10"/i~^ Mq, of which one-third is due 
to the backgroimd A— term. This is in the range 
of rich cluster mass estimates given by N. Bal- 
call ( in AQ). The agreement is surprisingly good 
since we have ignored halo-halo collisions, result- 
ing from the galaxies' appreciable velocities, and 
ignored density structure within the cluster. For 
N galaxies in a cluster, the total cluster mass is 
DC iV^/^ so that this estimate adequately represents 
poor clusters as well if we use e.g. N = 8. 

Finally, a similar argument must limit the range 



of applicability of the exterior Schwarzschild solu- 
tions in case I models; the total force on a test 
mass must vanish when r = Rq, the mean spac- 
ing between galaxies. We must introduce a A 
term. Then, for Rq > r > Ri, the force is / = 
-a'^a'^Ri (6r2)-i - Ar/3 where A = ^a'^a'^RiR^^. 
In effect the distinction between case I or Case 
II is whether or not the halo is limited by back- 
ground baryon density or field density (or by 
force-balance). 

In actual cases, such as treating the Local 
Group of galaxies, the convenience of assuming 
spherical symmetry should be replaced by the 
weak field approximations for the metric tensor, 
requiring continuity of the metric and its first 
derivatives across interfaces. In this c;ase, the 
equations of the interfaces between halos will be 
more complicated and A and A will be functions 
of position. A particular halo, such as that of the 
Milky Way may combine elements of both cases 
I & II. The rotation curve for the Milky Way 
should be different in the directions toward and 
away from Andromeda at large r because the ro- 
tation curve (determined by /) is A— dependent. 

The requirement A > is not dictated by the 
physics of the 0— Lagrangian but is required to 
meet boundry conditions we imposed by first sub- 
tracting off surrounding material. As used, it is 
a representation of Mach's principle since it ex- 
plicitly appears as a 'background' mean density 
induced by the proximity of other halos. 

3.2. Close Interactions 

Because the theory is nonlinear, wave interac- 
tions between close galactic halos are much more 
complicated than in simple linear potential the- 
ory. But, some simple features are easy to see. 
Suppose one considers a case in which one wants 
to represent (j) as the sum of two components, 
(j) = ri{r, t) + ^(r, t). For example, one may repre- 
sent the established LDM field of one galaxy and 
the other the influence of that of another passing 
galaxy The time of interaction between colliding 
galaxies is ^ 10^ yr, much longer than the char- 
acteristic scalar field vibration time ^ 10^ — 10^ 
yrs. The coupling between the fields induced by 
the cubic terms may be studied by writing equa- 
tion (1) in two parts: 
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4. Discussion 



[d\t - V2]r? = m-'rtil - 3^^) - mW\ (H) 
[d\t - V^]^ = m-'^il - iif) - m^e- (12) 

If 3^^, 3?7^ <C 1, the two components act as in- 
dependent s-wave solutions, or independent halos, 
not affecting each other except through gravita- 
tional interaction (see the Appendix). If the in- 
eqalities do not hold, then the mixed cubic terms 
cause ^ and rj to act like Mathieu functions, per- 
mitting resonant interactions to develop during 
the long collision time. We shall use a crude 
model. Suppose we can time average some terms, 
^^77 V'^S, ~^ ('7^)C) regarding these av- 

erages to be very slowly varying functions of po- 
sition. Then, e.g. a quasi-static solution for ^ 
is given by using w^jj = m?{l — 3(r?^)). Conse- 
quently, if (77^) > 1/3, a steady state solution for ^ 
is not allowed; the original dark halo of this galaxy 
will be modified and may not recover after the col- 
lision. One needs a close collision, centers passing 
within ^ 10 kpc, for this to occur. In the past this 
may have been more frequent; if so, one expects 
that an equilibrium population of dark halos will 
have (^^) < 1/3. Since mp/f may be considerably 
smaller than m when {rj'^) 1/3 it is possible that 
many of second class of galaxies, those with very 
long (convex) rise distances, Ri^, to be the result of 
recent halo interactions. [Another way of getting 
long rise distances is to require a — > 1, since m 
should be replaced by m for small r; see equation 
(2).] 

In such close halo-halo collisions, the struc- 
ture of the halos are also changed by the gravi- 
tational interactions between the halos themselves 
and with any ordinary matter trapped in the ha- 
los' gravitational wells. Section A. 5 gives esti- 
mates of the gravitational forces involved. It is 
possible that some dark halos have been stripped 
of much of their entrapped luminous material. 

We have avoided discussing t-waves, regarding 
them as infinitesimal. But they may not be, in 
regions where halos collide or in earUer epochs of 
the universe. Suppose ^, represent a collection of 
'strong' t-wave fluctuations; then, in their pres- 
ence, the stability of s-wave halos will be affected. 



In Section 2 we eliminated all more conven- 
tional potential sources for the LDM. This pa- 
per explores the idea that dark halo galaxies are 
fluctuations in a scalar field. The s-wave halo so- 
lutions studied are sucessful in being necessarily 
spherical in their inner regions and in producing 
the form of the observed gravitational accelera- 
tion. In its simplest form, the only parameter that 
can vary from galaxy to galaxy is the square of an 
amplitude, < 1, which controls the depth of 
the gravitational well of the LDM. Also, there are 
parameters & which have been determined 
from the Milky Way's rotation curve; we expect 
them to have these values for all galaxy halos at 
the present epoch 2 w if our Lagrangian is com- 
plete. 

Finally, there is a parameter Ri which influ- 
ences the gravitational attraction in the very outer 
regions of a halo. Its value may depend upon the 
density of neighboring galaxies. We have guessed 
that Ri >~ 300 kpc for an average halo, but 
a much smaller value might be needed in a few 
galaxies (Casertano & van Gorkom 1991) . 

4.0.1. Tests and Utility Of the Theory 

The most important observational tests revolve 
around the question: What physics determines the 
values of and (and of i?/)? If this is the infla- 
tion field then these parameters must have varied 
over large time scales and be functions of z. 

If is constant at the present cosmic epoch, 
then angular measurements of Rh in spiral galaxies 
would provide a new cosmic distance scale. At- 
tempts to establish such a distance scale would 
provide a test for the assumption is constant. 
One might be able to place limits on the variation 
of with time. If , e.g. oc (1 -|- z)^^^ then also 
oc (l + z)^^'^ and the rotation curves of galaxies 
at z « 2 would show larger amplitudes than those 
nearby. 

Because we may have omitted possibly impor- 
tant coupling terms to matter and radiation, it 
is possible that both and may have much 
different values in the immediate neighborhood of 
stars, compact clusters, or galactic centers. In 
principle, this possibility could be restricted by 
setting constraints upon the rotation of the lines 
of apsides in double star systems and in the Sag 
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A* -S2 system ( Schodcl et al 2002). Analyzing 
the velocity fields of colliding galaxies and com- 
pact galaxy clusters provides both a test and a 
utility of the simple theory. The non-linear terms 
should cause halo interaction forces stronger than 
the ones predicted by Newtonian theory and halos 
should feel a stronger than expected gravitational 
attraction to ordinary matter. 

Since baryonic matter can flow into the gravi- 
tational wells of LDM fluctuations and contribute 
additional source terms to the gravitational fields, 
the rotation curves of galaxies will be compos- 
ite. Use of this scalar theory to subtract out the 
LDM contribution should allow determinations of 
L/M for stars in the inner regions of galaxies. 
One test is therefore whether these determinations 
make sense. Suppose there are equal opportuni- 
ties for baryonic matter to have settled in most of 
the large LDM wells. Chose galaxies with similar 
colors for their central regions, so that {L/M) is 
the same for these galaxies' central region stars. 
Then the central luminosity within Rh = 5 kpc 
in field galaxies would be a function of a^, mea- 
sured by V)^. This would produce a Tully-Fisher 
relation but restricted to the very inner parts of 
galaxies, r < Rh- For low luminosity galaxies, the 
halo gravitational potential may actually be deter- 
minable in the center regions. Then, by Liouville's 
theorem, one would expect to observe star densi- 
ties oc exp [/? / / dr] . 

4-0.2. The Equivalent Cosmological Fluid 

In order to get a useful interpretation of the 
energy momentum tensor for (/>— field fluctuations 
regarded as individual dark matter halos, it was 
necessary to introduce reference energy and mo- 
mentum levels by subtracting from the basic (p 
field Tafc a energy-momentum tensor, r(A, A), 
where Kp — ^a^m?\ — A, np = ^a^m?X + A, 
W = m^/m, and f{X,A)f = 0. In crowded 
fields of galaxies, A represents an explicit under- 
lying background field density term used in de- 
termining the metric coefficient A. It seems this 
feature is unavoidable in any theory in which the 
source terms are derived from a Lagrangian for- 
mulation. Fluctuations in a field must be referred 
to a reference level. Because the reference energy 
density level in a Lagrangian normally is not in- 
cluded but must be specified in Einstein's equa- 
tion, one has a choice of either initially augment- 



ing the Lagrangian or of introducing T(A, A) as 
an additional required source term. In performing 
the large-scale averages of source terms required in 
cosmology, the intimate connection between the (f) 
field Tab and T{X, A) can be lost and two separate 
'independent' averages may appear. One wonders 
if in the remote past when fluctuations in the La- 
grangians for the weak and strong forces were more 
important contributors to the total energy density 
whether such A, A terms contributed significantly 
to the early cosmological A term. 

A fluid representation of the energy momen- 
tum tensor for a particular spherical s-wave so- 
lution when A = A = is not rewarding because 
Txx,Tyy,Tzz are not equal and vary from place 
to place; it is only by averaging over the surface 
of an entire sphere r = constant that a pressure 
can be defined: 3p = (T^^) + (Tyy) + (T,,) = 
{Ttt) = p- The same relation holds for averaged 
t- waves. However, for a cosmological fiuid, one 
consisting of a great many individual s-wave re- 
gions in a unit volume, one may not adopt the 
the same equation of state used for representing 
light, (p) = (/5)/3, because in this case, in general 
A, A 7^ 0. 1 am skeptical that spacial averaging can 
be dismissed as a trivial problem, for the link be- 
tween a source field and (some of) the cosmological 
constant A— term can easily be lost. The conser- 
vation laws (r"*" — T"'');(, = necessarily involve 
A & A, modifying the definition of the effective 
density and pressure. Therefore, the cosmic time 
evolution of the s-wave field density fiuctuations 
need not be oc 7?.(t)~^, where TZ is the cosmic scale 
factor, because the distribution of the halos them- 
selves must be considered. 

Suppose we assume we can replace galax- 
ies by halos, each of the same mass, M = 
{c^ /QG){a^)a^R, (see equation (9), where R is the 
effective maximum size of the representative s- 
wave halo. Suppose they were equidistant from 
one another with a number density -0. Assume 
that in the recent past, a was constant and halos 
were neither destroyed or created. If R were an 
assigned multiple of m~^, the mass of each halo 
decouples from the general cosmic expansion and 
the LDM contribution to the cosmic mass density 
would scale as Mtp oc TZ{t)~^ . For case II so- 
lutions,when the halo mass is oc i?o = V'~^^^/2, 
the LDM halo contribution and the A term would 
scale as TZ{t)~^. 
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The present niiniimim value for ij.) is probably 
a mass (or luminosity) weighted Schechter func- 
tion (see AQ, p.581), V = 0.014/1^ halos Mpc"^. 
Then, for case II, using Rq = /2, the val- 

ues of a^a^ found for the Milky Way galaxy and, 
M{Rq)iP+pa = i^LDMPc, one finds Qlom ~ 0.45, 
independent of h, with pA ^ 0.15pc- Since there 
can be halos not associated with luminous galax- 
ies, it is reasonable to consider ip for halos is ~ 1.5 
larger than that used; then the values of 17 would 
then be 4/3 larger. These are conservative esti- 
mates, for suppose we just guess that each galaxy 
has a halo mass ten times its luminoiis baryon 
mass, then ^ldm ~ ^^^baryon — 0.3. The possi- 
bility that Q.LDM > ^baryon+^other is appreciable 
and Pa from clustering LDM cannot be ignored. 
If so, the current cosmic expansion rate is then 
determined by the LDM with TZ{t) oc t. 

This estimate permits us to point out an in- 
teresting scenario. Suppose we take at present 
Mother ^ 1/3 (= {l/2)nLDM), Scaling in time as 
ordinary matter, cx TZ{t)~^ . Then going back in 
time to when 7^ = 1/3 one sees the reverse would 
have been true, Vlother ^ 2Qldm and then the 
scale factor would have had a different time depen- 
dence, TZ{t) oc t"^/^. This would mean that in the 
time between z = 2 and the present, the universe 
would have been observed to experience an accel- 
erated growth rate. Such an acceleration may have 
been observed. [For a discussion of the observa- 
tions of Perlmutter, Riess and their very many col- 
laborators see Perlmutter (2003).] The argument 
can be inverted: If these observations hold up, it 
would be reasonable to conclude Oot/ier oc TZ{t)~^. 

5. Summary of Theoretical Uncertainties 

What we have found is that the cj)— field vari- 
ations, as presented, can explain features of ob- 
served galaxy rotation curves. Coupling to other 
matter fields is not presently needed since a non- 
linear theory can produce its own source term; 
however, future observations could easily require 
inclusion of such terms. 

There are some computational investigations 

needed to verify that the time- varying s-wavc so- 
lutions will settle down quickly enough to be rep- 
resented by the steady state solutions. Also, the 
fluctuation spectrum in physics is represented by 
the action of a stochastic force on the RHS of 



equation (1). In order to specify the cosmology 
of the (/(—field, the time evolution of this force 
is needed. At the very least it has to be shown 
that the evolution of (j) variations from the ob- 
served spectrum of fluctuations seen in the cosmic 
microwave background can produce the presently 
observed galaxy clustering. In investigating this 
problem, one would effectively determine the time 
dependence of the (/>— field back to the time of pho- 
ton decoupling. 

Because of the uncertanties in knowing how 
& varied with time, and what role should 
be ascribed the t-waves in the very distant past, 
we simply do not know whether the scalar field 
proposed here is compatable with any of the many 
scalar theory suggestions for 'dark matter' or 'dark 
energy' which have already been proposed. Con- 
ventionally one sets for the spacially averaged 
scalar field for the dark matter a fluid represen- 
tation, p = {dt(j)f/2 - V{(j)) - ((V0)2)/6 and 
p = {dt(t>f/2 + v{(j)) - ((V</))2)/2, with the spa- 
tial derivative terms normally ignored. [ See Kolb 
& Turner (1990) and comprehensive reviews of 
dark matter theory by Peebles & Rata (2003) & 
Bernardeau (2003).] We require that the averaged 
spatial derivatives be included. For us, the cos- 
mological average for p is the same, (substituting 
{dr(t)f for (y(t)f and using V((/)) = -\m^(t>'^(l - 
02/2)), but p is different p = {dt(t>f/2 - !/(</)) + 
{{dr4'Y)/2, even when the necessary terms A & A 
are ignored. The two forms are not compat- 
able, even for t-waves. [ The differences arise 
from the fact we used in the Lagrangian a factor 
m°"ma rather than m^; one gets a different energy- 
momentum tensor from the conventional one when 
the variation SC/Sfjab is performed. We prefer our 
form of the Lagrangian because we do not get ba- 
sic changes in the algebraic form of the Lagrangian 
if we consider simple transformations in gat like 

Gab -dab-] 
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A. APPENDIX 



A.l. Representations of the s-Wave Solutions 

We consider the flat space field equation (1) with A = B = 1, requiring solutions be initially infinitesimal. 
There are two classes of solutions. Consider initially local mode excitations, (j) oc expi{u}t + k ■ r); then the 
dispersion relation for the modes become: 

= k^ -m^il- {(f>*(f>)). (Al) 

For > wc have t-wave solutions which remain infinitesimal. The second class, s-wave solutions, corre- 
spond to < 0; they are non-traveling waves and may grow to finite amplitude, their smallest values being 
determined by boundry conditions. The Fourier transforms (j){k, t) for the s-wave solutions must have the 
feature that (j)(k,t) ~ for fc^ > rn^ (since these Fourier modes always remain infinitesimal.) Consider the 
rest frame of a point treated as the origin. This restriction implies that <j) must fall off rapidly for r > 7r/2m 
and is effectively large only in the small central region. 

For simplicity wc consider ^ to be a real field. All s-wave solutions finite at the origin, can be written as 
a{t)+((){r,t). Then 

= -m2</.(l - P), (A2) 

where /(r, t)^ < 1 is a slowly varying bounded function [ This follows from taking the Fourier transform of 
V'^cj) and applying the mean value theorem.] All s-wave solutions must be of this form including the steady 
state solutions, p = (j)^. For solutions ^ as r — > oo, we require a = 0. The time dependence is then 
given by: 

dl<f> = mV(f - .^'). (A3) 

With suitable boundry conditions, some solutions are stationary and some oscillate^ with p representing a 

'mean' value of (fp' . 

We need consider only those solutions for which, at large t, the maximum values of (f^ decrease as r 
increases; then p must also decrease with r because it is bounded by the maxima and minima of cjp'. When 
the oscillations are of low amplitude, one finds / w (ps, one of the steady state solutions discussed below. 
More generally, we shall assume that in time both {p) and (fP' approach one of the steady state values 
(1)1. Since the frequency ^ rris/ P also decreases with r, the time variation of the solutions discussed below, 
is effectively confined to the central regions and the steady-state solution adequently represents the outer 
regions. 

A. 1.1. Radial Symmetry 

Only spherically symmetric solutions can dominate in an inner radial region, since the non-spherical modes 
cannot satisfy the restriction fctransijerse = ^(^ +l)/r'^ < until r is large; e.g. for m ~ 3 kpc, the t = 2 
mode cannot possibly be significant until r >'~ 7 kpc. If many other isolated halo regions of s-wave solutions 
exist they will cause nonradial interactions in each other's outer regions. Considering the case of all regions 
being of the same minimal size; then £ = 6 is appropriate for regions in contact and each would need a 
minimum core radius > 18 kpc for such interactions to be represented by s- waves; if this restriction is not 
true t-waves will be generated by the interaction. 

A.2. The Steady- State Halo Solution 

Equation( A2), with p = (fp, may be replaced by an integral equation using a Greens' function and 
regarding the nonlinear term as a source. In the general case, the formal steady state solution, can be 

^Thc boundry condition at t = is that (dt(f>)'^ is small, satisfying equation (A7) with < Oq < (/^). If this condition is not 
met, the solutions will exponentially decay. 
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represented by a convergent series expansion in the amplitude a < 1 at r = 0. One has (f>s = cl{4>o + a'^'^i + 
a'^(l)2 + •••)' where 0o = smmr/nir (plus angular terms.) for an isolated halo. The 0„ can be related to one 
another by e.g. 

Alternatively, a static (or time averaged) spherically symmetric real solution can be obtained by sucessive 
approximations. Put ^„ = {fn/r) sin ; then 



= y [m2(l - </>2_ J + a2,/„_i]i/2 dr, 



(A5) 



where = {a/m)drgn) is required as a constraint. To lowest order we use (po = (a/m.r) sin mr and df.,^f = 0. 
Succesive approximations amount to an expansion in the amplitude a at r = 0. Since it is probable that 
< 1/3, the next approximation should be sufficient: 

<^ « = (a/mr)[(l - a^){l - 0^)]"^/^ sin J m[l - 502/6]i/2 ~ {a/mr^{l - a^)) sinm J [1 - <Pl]^/^ dr. 

(A6) 

A. 3. Time dependent Halo Solutions 

There is a class of time dependent halo solutions which are related to the steady state solution, permitting 
us to estimate the time it takes the s- waves to develop. In equation (A2), for some time interval, replace 
by a time averaged value, (/^), to simplify the time dependence; then equation (A2) can be integrated: 

{dtcl^f = im2(/2)(02 _ al){aj - <p'), (A7) 

where ao(^) the value of (/)^(r) at the lower turning point and a^(r) = 2(/^) — aQ(r) is its value at the 
upper turning point. Using the notation a§ = (/^)(1 — h?{r)) for defining h^{r), the solution may be written 
as: 

02 = {p)H where H=[l + h^{r) sinO], (A8) 
with the time dependence given by the standard elliptic integral 

[2m''{f)]^{t-to{r))= f H-'/^de (A9) 

J-tt/2 

For small values of this is basically Kepler's equation relating the mean and eccentric anomalies. For the 
steady state solution, /i^ = 0. 

A solution still requires a value of ( j^). One may regard equations (2) & (7) as defining (Z^) in terms of 
(fy^, by sucessive substitutions. If this value is substituted into equation (A3), one can solve to get the actual 
form of (j)(r,t).^ An adequate approximation, for /i^ ~ 1, is to ignore the spacial variation of H and replace 
(f) by ^sV^j where now is the (approximate) steady state equation given in the form of equation (A6): 

(j)'^^H ■ {a/m) sin j m(l - a^cpiy/'^dr. (AlO) 



^The formal procedure is to use the specification of </>(r, t = 0) in equation (A2) to solve for /(r, 0); using this and the specification 
of dt4>{r,t = 0) in equation (A7), one calculates OQ^r). The integration, equation (A9), then gives the time evolution in the 
vicinity of t = 0. Because the time varying amplitude so obtained is r— dependent, (^{r, St) will have a different r— dependence 
than 0(r, 0); this causes f{r,St) to change, requiring the process to be iterated for the next time step. While / is bounded, 
it too will oscillate with time. Equation {A3) forces mode-mixing and some generation of t-waves which travel away from the 
halo. This produces some damping of the time dependent s-wave solution. The use of {P(r)) permits a reasonable estimate of 
the rise time and of the value of {(j>{r)'^), 
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For evaluating components of the energy momentum tensor we note time derivatives of s-wave solutions 
normally can be ignored since (jP' ~ m'^{p)4>^ < rn?(j)^ ~ {dr4>)^- 



= 4,/. - {B/Ar^)dr{r^dr(l)) - BVl(t> - m'^(i){l - (j)'^) + -[dt(t> dtHA/B) - dr<j) dr{B/A)], (A12) 

Vl^ = (r^ sinHy^lsmedeismedec^) + dl^<j)]. (A13) 



A.4. The Rise Time 

We now consider the rise time of these s-waves. From the rotation curves, we estimate m^^ ~ 3 kpc so 
that the time scale is measured in units of ~ 10'* yr. For the case Uq <C the rise time, ti, from ao 
to ai is then approximately, 

m/ti «ln[2V2/7ao]. (All) 

The initial variations in Oq oc Tu are assumed to be ^ 1 x lO"**, equal to the flux density variations seen 
in the cosmic microwave background. Our solutions are limited to r < Ri where (f>{Ri)^ ~ 10"**, for the 
solutions cannot be carried out to values less than the fluctuation levels from which they arose. This 
determines Ri « 100m~^^a^ ^ 170 kpc for = ^ and rise times ^ 1 — 5 x 10^ yr for the central region 
mr < 10. Therefore the characteristic times of oscillation and of growth to full amplitude are short compared 
to present galactic internal dynamic time scales, but comparable to the age of the universe at the time of 
photon decoupling. 

A. 5. Curvature and Restricted Two-Body Problems 

For a Schwarzschild metric, the wave equation is: 

where 

^'2(/^=V' »ni u) [sm 1/(^6) vsni lyugc/^; -r i^^^f, 

For the weak gravitational fields considered in this paper ,the coefficients A, B ^ 1 plus terms of order [Vh /c)^ 
and are static. The term dt<t> dt(B/A), for a time-dependent halo, produces damping with a characteristic 
inverse time ~ (V;i/c)^m and can be neglected. The term drcj) dr(B/A) is of order ^ {Vh/c)'^r^^dr(l)) and 
produces a very slight distortion of s-wave solutions centered at the origin. So, neglecting terms of order 
(Vh/c)^ the flat space wave equation is suflacient to evaluate contributions to the energy momentum tensor, 
in equations (4)- (6). 

However, the gradient term, oc dr{B/A), does determine the motion of the centroid of a distant small 'test' 
halo, one which makes a negligible contribution to the energy-momentum tensor. For suppose we change the 
independent variables in the field equation (t,x) (i,w), the local coodinatcs centered on the test halo, 
X = u(t) -|- w, and as before, set A, iJ ~ 1. In terms of the new independent variables, the wave equation 
becomes the normal one for flat space, df^ — V<p — m'^<p{l — cjp') + . . . = 0, providing one sets 

ld{B/A). . 

u=--^^|_r, (A14) 

Here u = |u|, two very small terms, of order u^rn^cf) and umcjcp have been neglected, and we have used a 

Taylor's expansion of rf^/"^'* assuming u ^ w. This is an equation of motion for the centroid of the distant 
test halo in the field of a central halo or star. [ The next order term in the expansion gives the tidal force 
on the test halo.] 

If the central object generating the metric is an ordinary baryonic mass or black hole, one has AB = 1 
and A~^ = —2GM/r; consequently the test halo experiences twice the gravitational attraction that a test 
mass would feel. If the central object is a s-wave halo, and w ^ m then A fa constant, since M oc r; 
the distant test halo experiences the same acceleration / a text mass experiences (see equation 76)). Since 
halo-halo interactions are wave phenomena; the halos do not really make separate contributions to the 
energy-momentum tensor. In close massive halo-halo interactions, the energy-momentum tensor will show 
interference phenomena and the wave equation, because of the non-linear terms,will also contain interference 
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terms. Consequently interactions between close massive halos will not be as simple as those involving 'test' 
halos. 
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